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ABSTRACT 

Context. Ap Lib is one of the rare low-synchrotron-peaked blazars detected so far at TeV energies. This type of source is not properly 
modelled by standard one-zone leptonic synchrotron self-Compton (SSC) emission scenarios. 

Aims. The aim of this paper is to study the relevance of additional components that should naturally occur in an SSC scenario for a 
better understanding of the emission mechanisms, especially at very high energies (VHE). 

Methods. We use simultaneous data from a multi-wavelength campaign of the Planck, Swift-UVOT, and Swift-XRT telescopes car¬ 
ried out in February 2010, as well as quasi-simultaneous data of WISE, Fermi, and H.E.S.S. taken in 2010. The multi-lambda emission 
of Ap Lib is modelled by a blob-in-jet SSC scenario including the contribution of the base of the VLBI-extended jet, the radiative 
bloh-jet interaction, the accretion disk, and its associated external photon field. 

Results. We show that signatures of a strong parsec-scale jet and of an accretion disk emission are present in the spectral energy 
distribution. We can link the observational VLBI jet features from MOJAVE to parameters expected for a VHE-emitting blob accel¬ 
erated near the jet base. The VHE emission appears to be dominated by the inverse-Compton effect of the blob relativistic electrons 
interacting with the jet synchrotron radiation. In this scenario, Ap Lib appears as an intermediate source between BL Lac objects and 
flat-spectrum radio quasars. Ap Lib could be a bright representative of a specific class of blazars, in which the parsec-scale jet lumi¬ 
nosity is no more negligible compared to the blob and contributes to the high-energy emission through inverse-Compton processes. 
Conclusions. 

Key words. Radiation mechanisms: non-thermal - Galaxies: active - Galaxies: jets - Radio continuum: galaxies - Gamma rays: 
galaxies - BL Lacertae objects: individual: Ap Lib 


1. Introduction 

Blazars are a distinctive class of active galactic nuclei (AGN) 
that are known to exhibit rapid variability, to lack bright emis¬ 
sion lines in their spectra, to show strong and variable polariza¬ 
tion at optical and radio wavelengths, and flat high-frequency 
radio spectra (Strittmatter et al. 1972). They are known to 
emit radiation across the entire electromagnetic spectrum. Their 
broad-band continuum energy distribution (SED), from radio to 
gamma-rays, consists of two distinct broad bumps: a low-energy 
component from radio through UV or X-rays, and a high-energy 
component from X-rays to gamma-rays. In most standard lep¬ 
tonic scenarios, the first bump is interpreted as the synchrotron 
emission of highly relativistic electrons in magnetized plas¬ 
mas, the second high-energy component as due to the inverse- 
Compton emission from electrons upscattering synchrotron pho¬ 
tons (synchrotron-self-Compton scenario, SSC) or photons from 
the ambient fields (external inverse-Compton, EIC). 

Blazars can be classified according to the frequency of the 
peak of the low-energy (synchrotron) SED component: 

- Low-synchrotron-peaked (ESP) blazars, consisting of flat- 
spectrum radio quasars (ESRQs) and low-frequency-peaked 
BL Lac objects (LBLs), have their synchrotron peak in the 
infrared regime, at < lO'^ Hz. 

- Intermediate-synchrotron-peaked (ISP) blazars, including 
LBLs and intermediate BL Lac objects (IBLs), have their 
low-energy peak at optical - near-UV frequencies at 
lO'^Hz < V, < 10'^ Hz. 


- High-synchrotron-peaked (HSP) blazars, which are almost 
all high-frequency-peaked BL Lac objects (HBL), have their 
synchrotron peak at X-ray energies with > 10'^ Hz. 

The majority of blazars (almost 80%) detected at VHE be¬ 
long to the HBL subclass. Their average SED as well as indi¬ 
vidual flares can in most cases be well described by simple SSC 
scenarios. However, a few cases of LBL and IBL sources, among 
which Ap Lib (PKS 1514-241), have been seen at VHE as well 
and raise specific questions. 

ApLib was first quoted as a variable star by Ashbrook (1942) 
before it was identified as the optical counterpart of the radio 
source PKS 1514-241 (Bond 1971; Biraud 1971). It had also 
previously been identified with a faint, compact galaxy (Bolton 
et al. 1964). Strittmatter et al. (1972) and Bond (1973) first sug¬ 
gested that Ap Lib was a member of the BL Lac class of ob¬ 
jects. To date, the best available measurement that demonstrates 
its very fast optical variability gives a variation rate of up to 
0.06 + 0.01 mag/hr in the V band (Carini et al. 1991). Optical 
polarization as high as 8.0% is observed (Stickel et al. 1994), 
but the average degree of polarization is quite low in the IR 
(Andruchow et al. 2005). 

Laint emission and absorption lines confined to the nucleus 
were detected, which led to a first estimate of a redshift of 
z=0.049 (Disney et al. 1974), which is confirmed by the final 
redshift release of the 6dL Galaxy Survey (Jones et al. 2009). A 
supermassive black hole mass of 10*^ was deduced from 

velocity dispersion measurements together with an elliptic host 
galaxy mass of by Woo et al. (2005). In the same 
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paper, Ap Lib, known as an RBL (radio-selected) source (Ciliegi 
et al. 1993), was named LBL for the first time. 

Ap Lib presents a rather flat radio spectrum from 8.4 to 90 
GHz with an average spectral slope of ?; = 0.76 {vFy oc ■/>) 
(Tomikoski et al. 1993). A 43 GHz polarization was measured 
by Lister et al. (1998), only one component located at ~ 1 mas 
to the core showed a significant polarization level of 7%, the 
core itself had a polarization level below 0.08%. In the same 
paper the radio jet morphology is described with a jet pointing 
toward P.A.= 171°, within 1 mas to the core, then turning to 
P.A.= 157°. More recently. Lister et al. (2013) reported direct 
VLBI measurements of jet flow kinematics and magnetic field 
properties. ' 

Ap Lib was first detected in X-rays by the Einstein X- 
ray Observatory (Schwartz & Ku 1983). The high-energy (HE) 
emission of the source (E > 100 MeV) was detected by differ¬ 
ent atmospheric balloons (Erye et al. 1971), but only in 1999 it 
was realized, thanks to EGRET on the CGRO satellite, that the 
gamma-ray source is associated to Ap Lib (Hartman et al. 1999). 
In 2010, the H.E.S.S. collaboration Hofmann (2010) reported 
very high energy (VHE) emission (E > 100 GeV) detection. 

Currently, only ten other LBLs or IBLs have been detected in 
the TeV energy range by lACTs, namely BL Lacertae, 3C66A, 
W Comae, S5 0716-H714, lES 1440H-122, PKS 1424-1-240, VER 
J0521H-211 and MAGIC J2001H-435. These sources are rather 
peculiar, for example BL Lacertae and VER J0521H-211 could 
have a different categorisation following their activity, IBL in 
low state and HBL in high state. The recent observation of an ex¬ 
tended X-ray jet by the Chandra satellite (Kaufmann et al. 2013) 
in addition to its VHE detection makes Ap Lib a unique LBL 
specimen. Moreover, pure SSC models have difficulties to repro¬ 
duce the properties of at least three of these sources, BL Lacertae 
(Ravasio et al. 2002), 3C66A (Joshi & Bottcher 2007), and W 
Comae (Acciari et al. 2009). Scenarios that include an additional 
gamma-ray emission by an external comptonization have been 
favoured for these sources. More recently, Raiteri et al. (2010) 
highlighted the existence of two SSC zones in the BL Lac SED. 

Previous attempts to model Ap Lib with simple one-zone 
SSC were not successful (Tavecchio et al. 2010; Sanchez et al. 
2012). The very broad high-energy bump of Ap Lib suggests 
that a more complex SSC model is needed to describe the 
emission scenario, as we show here. 

In Sect. 2 we describe the data we used and the host galaxy 
correction. In Sect. 3 we follow the method of Tavecchio et al. 
(1998) to test the relevance of simple SSC scenarios and to con¬ 
strain the B-6 domain of magnetic held and Doppler factor. We 
also show that it is difficult to fit the SED with these constraints 
with a pure SSC code. In Sect. 4 we further develop a multi- 
component SSC code based on the ”blob-in-jet” scenario by 
Katarzynski et al. (2001) by taking into account specific spec¬ 
tral components due to the parsec-scale jet and accretion disk, 
which were negligible in the cases of HBL sources, but could 
play a significant role in IBL/LBL objects. VLBI radio data are 
used in Sect. 5 to deduce additional constraints on the model 
that we propose and discuss in Sect. 6. We report in Sect. 7 that 
other sources present an SED somewhat similar to Ap Lib. They 
might form a specific subclass of blazars together with Ap Lib 
with a significant imprint of their VLBI jet on the SED even at 
non-radio frequencies. 


* http://WWW.physics.purdue.edu/astro/MOJAVE/ 
sourcepages/1514-241.shtml 


We name the various components of the jet as follows: 

- Extended jet: jet seen in radio and X-ray from pc to kpc scales. 

- Pc-jet (parsec-scale jet): base of the extended jet until 100 pc. 

- Inner jet: corresponds to the string of radio knots observed 
along the jet axis. 

Throughout the paper, we consider a cosmological constant 
Ho = 70 km.s *.Mpc ' and Qm = 0.3. 

2. Spectral energy distribution 

2.1. Data 

Since Ap Lib is variable, simultaneous or quasi-simultaneous 
data are collected to build the SED. We used the contempora¬ 
neous data recently published by Giommi et al. (2012), includ¬ 
ing data from Planck (from 17-02-2010 to 23-02-2010), Swift- 
UVOT (20-02-2010 and 22-02-2010), Swift-XRT (20-02-2010), 
and Eermi-LAT (January to March 2010). ApLib was observed 
by H.E.S.S. from 11-05-2010 to 10-07-2010 (Eortin et al. 2010), 
with a spectral slope of E = 2.45 + 0.20 (Sanchez et al. 2012). 
No strong variability being highlighted by Eermi observations, 
we consider that in spite of the temporal gap with the multi¬ 
wavelength campaign, these data can be used in the SED mod¬ 
elling. Non-simultaneous data were used to complete the gen¬ 
eral shape of the SED, but were not be taken into account in the 
modelling. We show non-simultaneous points only in the energy 
ranges uncovered by simultaneous data. We included archive ra¬ 
dio data from the NASA/IPAC Extragalactic Database (NED) ^ 
and Swift-BAT data from the ASI Science Data Center (ASDC) 
Multi-wavelength data are shown in Eig. 1. Throughout the 
paper, the host galaxy contribution was subtracted as described 
below. 

2.2. Contamination by the host gaiaxy 

Ap Lib is a nearby source, therefore its host galaxy has been 
resolved in the optical and near-infrared (e.g. Kotilainen et al. 
(1998); Urry et al. (2000); Pursimo et al. (2002); Hyvonen et al. 
(2007)). To model the SED, this contamination of the non- 
thermal emission has to be taken into account. 

Imaging studies of BL Lacs have shown that the host galax¬ 
ies have de Vaucouleurs type surface brightness profiles, as ellip¬ 
tical galaxies do. Moreover, R-H colours (Hyvonen et al. 2007) 
agree well with evolutionary population synthesis models, which 
predict for low-redshift, solar metallicity ellipticals R-H~2.4 for 
a typical age of 13 Gyr (e.g. Eioc & Rocca-Volmerange (1997); 
Bruzual & Chariot (2003)). To evaluate the host galaxy contri¬ 
bution to the observed radiation of Ap Lib, we thus used the el¬ 
liptical galaxy template spectra from PEGASE 2 (Eioc & Rocca- 
Volmerange 1999) to generate a synthetic spectram for a galaxy 
of age 13 Gyr and mass of (Woo et al. 2005), at 

the distance of the source. 

Then to estimate the host contribution in each UVOT and 
WISE Alter of the SED, we used the effective radius deduced 
from two-dimensional R-band photometric decomposition into 
nucleus and host galaxy, = 6.72” by Pursimo et al. (2002) 
and the aperture correction given in Eq. (4) of Young (1976) for 
a de Vaucouleurs profile of a spherical galaxy. 

In all subsequent figures the host contribution has been sub¬ 
tracted. 


^ http://ned.ipac.caltech.edu/ 
^ http://tools.asdc.asi.it/ 
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3. Constraints on a simple SSC scenario 

The unusual width of the high-energy bump in the SED makes 
modelling with pure SSC scenarios very challenging. No SSC 
model published on this source does account for VHE data 
(Tavecchio et al. 1998; Sanchez et al. 2012), and other radiating 
components are probably needed to describe the Ap Lib spec¬ 
trum. However, before considering more complex models, we 
first analyse the limits of pure SSC scenarios for AP Librae in 
this section. 

Under the SSC assumption, we can derive physical con¬ 
straints on the source from the SED shape and the variability 
properties and border the B-5 domain for the magnetic field B 
and Doppler factor 6. 

We considered a simple stationary homogeneous one-zone 
SSC model (Katarzyhski et al. 2001). The SSC model describes 
the emitting region as a spherical plasma blob of radius R, 
filled with a tangled magnetic field B, that propagates with a 
bulk Doppler factor 6 inside the jet. The low-frequency and 
VHE bumps are interpreted as the synchrotron and IC emis¬ 
sion from a population of relativistic leptons (e+) assigned to 
a single-electron population for simplicity. The primary elec¬ 
tron energy distribution, between Lorentz factors 7 ^,,, and ymax, 
is parametrized by a broken power-law function, with normal¬ 
ization factor K (defined as the number density of electrons at 
y — 1 , in units of cm“^) and indices n\ below and n 2 above a 
break energy yhreak- The absorption by the extragalactic back¬ 
ground light is taken into account by adopting the model of 
Eranceschini et al. (2008). This type of SSC model has typically 
eight main free parameters and thus requires high-quality data 
over a wide spectral range to be well constrained. 

The size of the gamma-ray emitting zone R can be derived 
from the minimum variability time scale f,,ar following the ex¬ 
pression R - tyar^cjil + z). To our knowledge, the fastest re¬ 
ported variability of the source is tyar = 1-3 days. A luminosity 
variation of a factor about 1.8 was observed in ~ 24h in the V 
band between 1989 March 16 and March 17, and was extrapo¬ 
lated to a doubling time scale (Carini et al. 1991). 

The shape of the observed SED can be used to derive ob¬ 
servable quantities by imposing further constraints on the phys¬ 
ical parameters (see e.g. Tavecchio et al. (1998)). The frequency 
and luminosity of the synchrotron (v^, L^) and Compton (v^., Lc) 
peaks are estimated by interpolating the low- and high-energy 
bumps of the SED by a cubic polynomial function (see Eig. 
1). Because the low-energy component is highly constrained by 
the data, the synchrotron peak is well determined. Although the 
Compton component peak is poorly constrained by the avail¬ 
able data, it can be considered to lie between Swift-Bat and 
H.E.S.S. data. Below and above the synchrotron peak, the spec¬ 
trum is quite smooth and can be approximated with a power law 
{Fy oc v") with indices ai = 0.29 and a 2 = 1.58, respectively. 
Adding the highest observed frequency in VHE, Vy, Table 1 sum¬ 
marizes the eight observational constraints provided by the data. 

Erom the condition of transparency of gamma-rays to pair- 
production absorption, a strong constraint on the minimum 
Doppler factor is obtained, 5 > 10.3. For the values of the 
peak frequencies and slopes in Table 1, the KN limit implies 
^KN < 0.29, which is lower than the transparency limit. We 
can thus consider that the emission scenario is dominated by 
the Thomson regime. The synchrotron and Compton peak ratio 
gives y^. The ratio of the total luminosity of the synchrotron peak 
to that of the inverse Compton is directly related to the ratio be¬ 
tween the radiation and magnetic field energy density inside the 
emission zone. This allows deriving a relation between the mag- 



Fig. 1. Interpolation of the Ap Lib SED by two-third de¬ 
gree polynomial functions, showing the synchrotron compo¬ 
nent (red) and the inverse-Compton component (blue). Black 
data: Simultaneous and quasi-simultaneous data, dots: Planck, 
Squares: WISE, stars: Swift-UVOT, diamonds: Swift-XRT, 
crosses: Fermi-LAT, grey dots: Archive data for informational 
purposes only, grey bow-tie: H.E.S.S.. The red rectangle is 
the uncertainty area of the Compton peak position. Only quasi- 
simultaneous and H.E.S.S. data are used for the interpolation. 


Table 1. Observables deduced from the SED shape. 


Symbol 

Value 

Unit 

Ls 

1.25 X 10'*'* 

erg.s^*. 

4 

[5.1 X 10''^L3x 10““] 

erg.s^' 

Vs 

5.9 X 10'^ 

Hz 

Vc 

[2.0 X 10‘^5.0x 10^'*] 

Hz 


4.0 X lO^® 

Hz 

a, 

0.29 


ai 

1.58 



netic field and the Doppler factor, and with estimating the source 
dimension, a domain of parameters such as 10.3 < 5 < 100 and 
72.5/(5^ < B < 344/(5^ in Gauss. 

Figure 2 illustrates the typical best fit that can be obtained by 
respecting the parameter domain B-d. The model parameter val¬ 
ues are given in Table 2. The shape of the SED imposes the need 
for a steep synchrotron slope 112 above the peak, so that the X- 
rays around 10'^ Hz are not overproduced. This constraint alone 
prevents the Compton component from reaching the observed 
VHE flux, regardless of the choice of the other free parameters, 
even considering second-order scattering, where incident pho¬ 
tons are the first-order inverse-Compton photon density. 

In addition, other SED features are not explained by the one- 
zone SSC scenario, such as the low-frequency radio emission, 
the flatness of the Fermi spectrum, and the change of slope in 
radio-mm at 250 GHz. The reducedon a simple power law fit¬ 
ted on radio data confirms the significance of a change of slope in 
radio-mm range - 3-7). The synchrotron bump shape 

does not account for the far-UV data, there is a clear excess 
above the model. This excess is usually associated with the ther¬ 
mal emission from an accretion disk. 

We demonstrate in the following section by taking into ac¬ 
count several physical effects that are neglected in basic SSC 
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log(E [eV]) 

-6 -4 -2 0 2 4 6 8 10 12 14 



Fig. 2. Best tentative SSC modelling of the SED, within the B-d 
constraint. Blue lines: synchrotron and SSC emission from the 
blob, blue dashed line: second-order SSC emission. Parameter 
values are given in Table 2. Basic SSC scenarios are unable to 
reproduce the VHE fluxes (see text). 


Table 2. Values of physical parameters used to model the SED 
shown in Eig. 2, consistent with the constraints deduced from the 
observables given in Table 1. 0 is the angle of the blob direction 
with the line of sight, ni and n 2 are the first and second slope of 
the electron spectrum, R is the radius of the emitting region. 


Parameter 

Value 

Unit 

6 

12 

- 

e 

1.0 

deg 

K 

1.5 x 10“* 

cm“’ 

«i 

2.0 

- 

«2 

4.5 

- 

7mm 

30 

- 

ymax 

1.0 x 10’ 

- 

ybreak 

5.5 x 10’ 

- 

B 

0.1 

G 

R 

1.8 x lO**^ 

cm 


scenarios, but are expected to naturally occur in sources, that all 
of these SED features can be integrated in a multi-component 
self-consistent scenario that is able to reproduce the SED up to 
very high energies. 


4. Multi-component scenario with jet and disk 

After rejecting the simple one-zone SSC scenario, we focus here 
on the relevance of other components that can play a significant 
role in the emission processes. Up to now, the influence of the 
extended radio jet has been mostly neglected when studying the 
multi-lambda emission of VHE blazars, except in a few stud¬ 
ies where it is invoked to explain the very low frequency radio 
spectrum. It is commonly assumed that its contribution to the 
total radiative output is negligible compared to that of the VHE 
blob, which is historically justified by the HBLs studies in which 
the SSC emission from the blob dominates the SED. However, 


X-ray observations of a bright extended jet by Kaufmann et al. 
(2013) suggest that this is not the case for the blazar Ap Lib. The 
two non-thermal components in the radio-mm spectrum deduced 
from the change of slope in the SED at 250 GHz (see Sect. 3) 
supports this view. One component is allocated to the blob and 
another to the jet. They are associated to SED peaks in the ranges 
[8.82 X 10\ 5.55 x 10^] GHz and [1.37 x 10^ 1.36 x 10"^] GHz. 

It is relevant to consider two populations of relativistic 
electrons in the non-thermal emitting zone of Ap Lib, which are 
related to the jet and to a blob embedded in it and contribute 
to the SED. In this case it becomes necessary to take into 
account the radiative components related to the jet particle 
population, mainly the jet SSC contribution and the inverse 
Compton of the blob particles on the jet synchrotron radiation. 
Thanks to radio VLBI data, which allow a detailed kinematic 
description of the jet as discussed in Sect. 5, this new compo¬ 
nent provides some constraints on the jet parameters without 
additional degrees of freedom in the global picture of the source. 

We propose a revisited blob-in-jet model, taking into account 
an inhomogeneous jet describing the non-thermal radiation from 
one peak seen in the synchrotron bump, including SSC and mu¬ 
tual inverse-Compton effects. The new model implemented here 
is based on the work done by Katarzyhski et al. (2001) and 
Katarzyhski et al. (2003), which was applied to HBL sources 
where effects related to the pc-jet were negligible. Hereafter we 
describe only the further developments and improvements in the 
scope of our study. We took into account the absorption of the 
VHE blob emission by pair production due to photon-photon in¬ 
teraction with the jet radiation, and the Comptonization of the 
jet radiation by blob particles. The blob moves with a Doppler 
factor higher than that of the jet, which induces Lorentz trans¬ 
formations between the blob and jet frames. To complete the 
picture, we also include a basic description of the disk emission 
and scattering of its light by the broad-line region (BLR), which 
probably plays a role in the SED of LBL sources such as Ap Lib. 


4.1. Basic jet modelling 

The jet is represented by a cone discretized in cylindrical slices 
(see Eig. 3). Absorption and emission coefficients are calculated 
for each slice. The radiation transfer is computed in the direction 
of the jet propagation. We express the intensity for each slice as 

Iijet(y) = [l “ eXp(-T;,;ef(v))] , (1) 


where TijJv) = LijethjeAv) is the optical thickness, Jijetiv) 
is the electron self-synchrotron emission coefficient, kijetiv) the 
absorption coefficient , and v the frequency in the jet frame. 
Eurther absorption by pair production was also considered as de¬ 
scribed hereafter by Eq. 6. 

The electron density Njet evolves along the jet with the slice 
radius as 

N,je, ^ (2) 

According to our conical description of the jet, we can de¬ 
scribe the expansion speed. 


diRi.jetj 

dt 


^ 13 jet tana = /3 jet 


DiJet 


(3) 


with (3jet the jet velocity assumed constant, and Dtjet the dis¬ 
tance of slice i to the central engine. We consider that the pc-jet 


4 











O. Hervet, C. Boisson and H. Sol: Linking radio and gamma ray emission in Ap Librae 


i = 50 



Fig. 3. Geometry of the jet. The jet is divided into 50 slices 
whose radius Rjet and width Lje, evolve along the jet accord¬ 
ing to the opening angle a. The shades of blue show the particle 
density and magnetic field decrease along the jet. 


4.3. External jet radiation on the blob 

The high-energy particle density in the blob is much higher than 
that in the jet since the so-called blob is by definition a compact 
zone of efficient particle acceleration processes due to shocks, 
turbulence, or magnetic reconnection, with a large population 
of relativistic electrons. The inverse-Compton emission of the 
blob particles on the local jet radiation can be significant, while 
the inverse-Compton emission of the jet particles on the blob 
radiation remains negligible. To describe this effect, we calcu¬ 
lated two radiation fields from the jet, one behind the blob and 
one in front of it, calculated via a radiation transfer in the direc¬ 
tion of (resp. opposite to) the jet propagation. The sum of these 
two fields acts like an almost isotropic field in the blob region, 
and the corresponding jet frequencies and intensities in the blob 
frame can be simplified as follows: 


expansion is free with an expansion velocity d(Rijet)ldt of the 
order of the Alfven speed. 

Then we can express the magnetic field of the slice i: 


Following Eq. 2, we can write 


R. 


‘Jet 


D 


ijet 


^i.jet — ^IJet 


D 


I Jet 


Di. 


jet 


(4) 

(5) 


This corresponds to the variation of the magnetic field along 
the jet as actually observed for a sample of VLBl radio jets by 
O’Sullivan & Gabuzda (2010). 


4.2. Blob radiation absorption by the jet 


We assume that Rbiob ^ Rjet, such that the jet radiation ab¬ 
sorption by the blob is considered negligible. However, the blob 
radiation absorption by the jet must be taken into account. In 
the high-energy range, the blob and jet photons can interact and 
create e^e~ pairs. For each slice we used the approximation of 
Coppi & Blandford (1990) in cylindrical geometry to estimate 
the pair absorption. We define the optical thickness of a jet slice 

Ti,yy(€c) by 

1 3 

r'i.yyi.ecj — 0.2.Crf riijl /6cj~LiJet- (6) 


In the source frame Vbhb > Vjet, and according to the speed 
transformation law, we can express the blob velocity in the jet 
frame by 


yb-Vjet 
" , ybVje, ’ 

and thus express Fj^ and 6'^ as a function of V^. 


(7) 


In the jet frame, slices in front of the blob receive a 
blueshifted radiation from it. We can write the blob frequen¬ 
cies and intensities transformations along the line of sight as 
x'ij - Vbd'i^ and = IbS'j^ and handle the absorption of this radi¬ 
ation slice by slice along the jet with frequencies and intensities 
of the blob expressed in the observer frame as 


Vobs = ^*^(1 +J> 

(8) 

u. = /;(!) (i+z)-'. 

(9) 


^'jet = ^'jet'^jet (10) 

~ ^ j" ‘^^^Jet^Jet ~ ^jet^jet- (H) 

This component of external inverse-Compton on the jet radi¬ 
ation due to a blob-jet interaction is found to play an important 
role for the VHE fluxes of Ap Lib (see Fig. 6). 

4.4. Disk and BLR radiation 

As we have seen in Sect. 3, the UVOT spectrum suggests an 
additional radiating UV component. In AGNs, such a compo¬ 
nent can be most naturally attributed to the thermal emission 
of the accretion disk. The shape of the SED shows that the 
UV flux density is at least higher by a factor 10 than that of 
X-rays, therefore we estimate that the disk luminosity is widely 
dominated by the thermal bump that we approximated here by a 
black-body at a given temperature. 

At higher energies the Fermi spectrum is rather flat. This 
feature cannot only be described by the blob and jet inverse- 
Compton components that are due to the strong constraint on the 
steep UV synchrotron slope (see Sect. 3). It suggests that an ad¬ 
ditional component radiates in the Fermi range that we can nat¬ 
urally identify as an additional external inverse-Compton emis¬ 
sion. 

As a result of the high blob Lorentz factor, the disk radiation 
interaction with particles in the blob frame is too highly red- 
shifted to be significant, therefore we expect that this inverse- 
Compton emission instead originates from the BLR disk that 
reprocesses radiation on the blob. To model this radiation, we 
approximated the accretion disk to a point-like source {Dbiob » 
Rdisk)- We assumed that the blob is located inside the BLR and 
propagates into an external isotropic photon field that is initially 
emitted by the disk, but is reprocessed by the BLR. The intensity 
of this field can be expressed as 


hlRiy, Tdisk) - TblR 


Rdisk Lplanckiy, Rdiskj 

AttRIlr (o-l7T)T^i,k 


( 12 ) 


where cr is the Stephan-Boltzmann constant, Ipimck is the 
Planck intensity, and tbir is a function of the reprocessing ef¬ 
ficiency and the covering factor. Lorentz transformation of this 
radiation field applies in the same way as in Eqs. (10) and (11). 


5 





































O. Hervet, C. Boisson and H. Sol: Linking radio and gamma ray emission in Ap Librae 


5. Constraints on blob and jet from VLBI radio data 

We have shown the influence of the blob and jet on the radiative 
output of Ap Lib. Now it is essential to constrain their physical 
parameters as much as possible so as not to increase the number 
of free parameters in the global MWL scenario. We rely for this 
on Lister et al. (2013) and focus especially on the data obtained 
by the MOJAVE programme at 15.4 GHz on December 26,2009, 
the closest epoch to the MWL observation campaign of February 
2010 . 


Table 5 of Lister et al. (2013), three of these knots show a 
clear superluminal motion with mean apparent velocities j3app of 
6.03(+0.34)c, 6.41(+0.19)c, and 6.04(+0.43)c. In the following 
we use the average value Papp = 6.16 c. We can link the apparent 
velocity to the real velocity p with 




I sin 0 
-h cos 

\Papp 



(15) 


where 6 is the angle between the jet direction and the line of 
sight. 


5.1. Radio core features 

The Ap Lib radio core on December 26, 2009 shows an ellipse 
shape with the major axis oriented at 341° on the sky plane as 
described in the Table 4 from Lister et al. (2013) . The core flux 
density detected is vFvcore - 1-52 ± 0.08 x 10“'^ erg.s“'.cm“^ 
(see Fig. 6 ). 

This core emission is strongly associated with the standard 
extended emission zone in the radio range of the SED. In our 
model we associated this radio core emission with the base of the 
simulated pc-jet that strongly emits in radio. Assuming a conical 
jet description and small angle approximation, we can thus con¬ 
strain the modelling parameters of the pc-jet via the radio core 
features. We assume that the half-minor axis of the ellipse cor¬ 
responds to the maximal radius of the pc-jet base, which has the 
core flux density, namely Rcore = 1.4 + 0.6 x lO'^cm {Dp = 218 
Mpc). This means that we have to resolve the following equation 
system: 


Rfi.jet — Rcore^ ^ — f'^Fv)core^ (13) 

i=l 

where YIi=\ ths sum over the n pc-jet slices contributing to the 
core emission (see Fig. 3). 

The solution of this equation system is strongly model 
dependent. To And the appropriate n slice, we adjusted the 
pc-jet parameters taking into account the modelling constraints 
given by the SED shape. With these constraints the pc-jet space 
parameter is strongly reduced. When we fulfil these constraints, 
we can estimate the half-opening angle of the pc-jet a - 0.4° 
and the size of the core emitting region Li,jet = 9 pc for 

a radius= L9xl0*^cm consistent with the measured value. 

Moreover, we assume that the major axis of the radio core 
ellipse Mucore - 8.6 + 3.4 x 10'^ cm corresponds to the size of 
this emission zone projected on the sky plane. With a conical jet 
seen at small angle, the length of the core projected is L^ore - 
Mucore - Rcore - 6.7 + 3.4 X lO'^cm, and so the viewing angle ft 

6 = arctan ( ] = 1.4° + 0.7°. (14) 

\ 2ji=1 Li, jet I 

Because of the adiabatic expansion of the pc-jet, the 9 pc 
core widely dominates the global emission of the total simulated 
pc-jet. The pc-jet emission between 9 and 100 pc contributes 
significantly to the global SED only in the radio range. We de¬ 
cided to fix the pc-jet length at 100 pc because the emission of 
the larger scale jet is negligible in the SED. 

5.2. Radio knot veiocities 

During the observation of December 26, 2009, seven radio knots 
have been identified in addition to the core emission. In the 


The Doppler factor 6, 


I - P cos 6 ’ 


(16) 


has a mean value of d = 21 . 8+4 4 , which is consistent with 
typical VHE blazar Doppler factors. Within this scheme, the ra¬ 
dio knot velocity is similar to that of the simulated blob in the 
blob-in-jet VHE scenario. 


5.3. Position and size of radio knots 



knots - core distance [pc] 

Fig. 4. Core distances and radius of radio knots observed by 
MOJAVE for an angle with the line of sight 6 of 1.38°. Grey 
dots show the referenced knots from August 18, 1997 to March 
5, 2011. Blue dots show radio knots of December 26, 2009. The 
black line is a linear regression that we use to characterize the 
knots expansion angle (j). The red dashed line at 100 pc marks 
the length of the simulated jet. 


To determine the distance to the core and the size of ra¬ 
dio knots, we considered all VLBI knots observed by MOJAVE 
from August 18, 1997 to March 5, 2011 referenced in the Table 
4 of Lister et al. (2013). The knot-core distance Dk was de¬ 
duced from the apparent (projected) knot-core distance Dhproj 
as Dii = Dk^projl sin 6 , and assuming spherical knots and small 
observation angle such that R^ - Rk.proj- The evolution of the 
radius with the core distance can be fitted by a simple law 
of proportionality R^ cx with a high level of confidence 
(X^Idof - 0.98). We can then estimate the half-opening angle of 
the knot expansion (p from a linear regression in the knot radius 
as a function of their distance: 

(^ = arctan = 0.10°. (17) 
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Fig. 5. Sketch showing the knot expansion angle cj), the opening 
jet angle a, and the angle with the line of sight 6. Not to scale. 


This yields (p < a < 9, and the jet does not strictly point in 
the observer direction. Because the observation angle 9 is very 
close to the jet aperture angle a {9 - a < 1°), we assume that 
radiation transfer formulae are still valid for such a small mis¬ 
alignment. 

This geometry describes two jets with two different open¬ 
ing angles, one innermost jet corresponding to the knot expan¬ 
sion and an external jet with a wider opening angle (see Fig. 5). 
This interpretation is consistent with the observations of Perucho 
et al. (2012) and Lister et al. (2013), who deduced a thin ribbon¬ 
like structure embedded within a broader conical outflow in sev¬ 
eral VLBI quasars. 


6. SED modelling and discussions 

From the modelling of the different components presented in 
Sect. 4 and the constraints on physical parameters deduced in 
Sect. 5, we can now generate the SED of Ap Lib in Fig. 6 with 
the parameters of Table 3. As we can see in Fig. 6, the unusual 
position of the jet and blob synchrotron peaks provide a very 
good fit of the lower part of the SED, in addition to a good re¬ 
production of the HE and VHE spectra. 

Erom modelling the UVOT data with a black body, we de¬ 
duce a temperature of 3.2 x 10^ K and a luminosity of 5.0 x 10"^^ 
erg.s \ which match the usual features of an AGN accretion 
disk. 

The flux density of each VLBI radio knots (blue dots in Eigs. 
4 and 6) is well reproduced by the synchrotron emission of the 
blob unabsorbed by the jet, and the core flux density (red dot) by 
that of the pc-jet. This gives us an overall view of the emission 
processes. The radio core emission comes from the base of the 
jet. We cannot see the radio emission of the emerging blob be¬ 
cause of the thickness of the core. But at a large distance from the 
core, the jet becomes optically thin in radio, and moving blobs 
can be seen as radio knots. Thus, the SED modelling supports 
the idea that the simulated VHE blob will become a radio knot. 
According to this model, the jet becomes radio transparent at 
15.4 GHz at about 6.5 pc from the core. This distance projected 
on the sky plane corresponds to a gap fro, the core of 0.24 mas. 
The average separation of the closest radio knots from the core 
reported by Lister et al. (2013) is roughly 0.85 mas, which is 
consistent with our proposal. In the following part of this section 
we examine different aspects for characterizing the link between 
the VHE blob and the radio knots. 


Table 3. Physical parameters used for the SED modelling pre¬ 
sented in Eig. 6. L is the size of the jet, Dhiob-BH is the distance 
between the blob and the central engine. 


Blob parameters 

Value 

Unit 

6b 

22 

- 

% 

1.4 

deg 

Kb 

2.0 x 10= 

cm-2 

ni 

2.0 

- 

111 

3.6 

- 

'Yminlb 

600 

- 

'Yma.xlb 

4.0 x 10* 

- 

'Ybreak/b 

8.0 x 10^ 

- 

Bb 

6.5 x 10-2 

G 

Rb 

6.2 x 10‘* 

cm 

Jet parameters 

6 jet 

10.0 

- 


50.0 

cm-2 

^jet 

2.0 

- 

ymax! jet 

1.65 x lO'* 

- 

Bl.jet 

8.0 x 10-2 

G 

^ 1 Jet 

2.5 x 10‘* 

cm 

L 

100 

pc 

a 

0.4 

deg 

^blob-BH 

7.9 x 10‘® 

cm 

slices 

50 

- 

Nucleus parameters 

Tdisk 

3.2 x 10^ 

K 

^disk 

5.0 x 10^2 

erg.s-' 

Rblr 

7.9 x 10‘* 

cm 

T BLR 

3.5 x 10-2 

- 


6.1. Cooling time 

The characteristic synchrotron cooling time Tcooiij) is 


Tcooliy) = 


3meC 

Ao-jUg fP-y’ 


(18) 


with cTj- the Thomson cross-section and U'g - the 

magnetic energy density. 


We obtain a characteristic synchrotron cooling time of 12.8 
hours for the very high energy of the blob (ymax - 4.0 x 10®). 
This time is much shorter than the year measured between the 
various radio knot ejections. In a purely ballistic scenario, where 
a high-energy blob is emitted with an injection time of ti„j - 
Rbic ^ 57h, then moves along the jet and emerges as a radio 
knot, we should detect very high energy flares only during the 
rare ejections of knots. However, the source appears stable at 
very high energy over years, which indicates that the particle 
acceleration must be continuous. 

The 15.4 GHz detection frequency of the radio knots is very 
close to the minimal blob frequency (see Eig. 6), so we assume 
that the mean Lorentz factor jr of the particles emitting in radio 
corresponds to the blob y„,i„ (see Table 3). The associated syn¬ 
chrotron cooling time is Tcooiijr) - 9.7 years. 

However, the knots with apparent movement detected by 
MOJAVE occur between 209 pc and 470 pc to the core on 
December 26, 2009 (see Eig. 4), which means assuming a con¬ 
stant velocity that they were ejected about 700 to 1500 years ago, 
a time much longer than the cooling time. Thus the particle ac¬ 
celeration should occur not only at the base of the jet, but also 
constantly at large distances from the core. 
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log(E [eV]) 

-6 -4 -2 0 2 4 6 8 10 12 14 



Fig. 6. Multi-component SSC modelling. Back data: simultaneous or quasi-simultaneous observations, grey data: non-simultaneous 
observations, blue dots: VLBI radio knots, red dot: VLBI radio core. Values of physical parameters are given in Table 3. The EBL 
absorption based on the model of Franceschini et al. (2008) is taken into account in this modelling. Because of the low redshift, the 
EBL unabsorbed and absorbed spectra are very close, we do not show the unabsorbed spectra for sake of visibility on the SED. 


The poor angular resolution at very high energy does not al¬ 
low determining where the gamma radiation zone lies. Two dif¬ 
ferent scenarios can be considered: 

(1) Either the gamma emission zone is stationary at the base 
of the inner jet and radio knots are created during large-scale 
perturbations due, for example, to a shearing between inner and 
extended jet. 

(2) Or the gamma emission zone propagates through the jet 
at the radio knot velocity and continues to radiate at high energy 
above the pc scale. 

The radio cooling time estimate, which requires an efficient 
acceleration mechanism at large distances, tends to favour the 
second scenario. Information on the variability pattern at high 
energy on long time scales is needed to disctinguish between the 
two options. The detection of a gamma flux periodicity corre¬ 
sponding to the radio knot periodicity (~ year) would support 
the second one, but this is beyond the scope of the present study. 

In the following part we study the jet energies in detail to 
characterize this continuity between the radio knots and the VHE 
blob. 


6.2. Cold particle population 

Figure 4 suggests a linear expansion of radio knots with their 
propagation along the jet. Assuming that the high-energy emerg¬ 
ing blob follows the same expansion as the observed radio knots, 
we can approximate its expansion velocity with the Alfven ve¬ 
locity VA,b as 

VA.fe tan 0 ^ 1.74 X lO^^c. (19) 

The Alfven velocity depends on the mass density p and the 
magnetic field B, we can deduce the blob mass density pt, from 
Table 3: 

Ph = ^ (—) - 1.2 X 10^*®g.cm^^ (20) 

47r \vA,bj 

That the proton mass is about a thousand times heavier than 
the electron mass implies that the blob mass density is widely 
dominated by protons. We can deduce the proton density Nb.p, 
assuming that the environment is electrically neutral, and ne¬ 
glecting any significant cold population composed of e /e^ pairs, 
which would annihilate very quickly, 

Nb„^ — ^l.Ax lO^cm^l (21) 

nip 
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We determine the density of non-thermal electrons Nb,e,NT 
by integrating over the non-thermal energy distribution spectrum 
Neiy): 

Nh.e,NT ^ f A^i,£(r)'^r - 1-8 X lO^cm^^, (22) 
Jymm 

with values of and ymax presented in Table 3. The cold 
electron population for this neutral environment is 

Nh,e,coid = Nb,p - Nb,NT = 7.4 X lO^cm-^ (23) 

In this way, we characterize the three blob components, 
namely the non-thermal electrons, the cold electrons, and 
the cold protons. The cold population of electrons appears to 
represent 99.8% of the total electron population. 

The same calculation for the pc-jet at its base provides an 
Alfen velocity vaj - 6.9 x lO^^c corresponding to the mass 
density pj - 7.7 x lO^^'g.cm^^. The proton density is Nj^p - 
4.6 X 10^cm“^, very close to the cold electron density Nj^e,coid — 
4.5 X 10^cm“^ , which represents 97.2% of the total electron 
population Nj^e,tot- The calculated non-thermal electron density 
2iNj^e,NT = 1.3 X lO^cm^^. 

In both the blob and the pc-jet, the cold particle density 
widely dominates the density of non-thermal particles, which 
does not contradict our assumption of continuous evolution be¬ 
tween the VHE blob and the radio knots. The wider aperture of 
the pc-jet results in a cold particle density lower than that of the 
blob. 

6.3. Energy budget 

Now we can examine energetic questions. Following the descrip¬ 
tion of Celotti & Ghisellini (2008), we can determine the power 
of various components through a section nR^ of the jet by the 
equation 

Pi - nR^T^dU'^), (24) 

where U'. is the energy density of the ith component in the 
comoving frame. The different components of power represent 
the magnetic field (Pb), the non-thermal electron population 
iPe.NT), the cold electron population {Pe.coid), the cold proton 
population {Pp), and the radiation (Pr) contribution. 

The associated energy densities are 

U'^ (see Eq. 18), 

U'e,NT = f yNe{y)dy, (25) 

where Ne(y) is the non-thermal energy spectrum, 

K.coid = rneC^Ne,coid and U'^ = mpC^Np, (26) 

with the cold densities Ng^coid and Np calculated in Sect. 6.2, 

t/; = — /;,ufv', (27) 

c 

where /', is the surface intensity. 

We calculated the powers of the blob and of the pc-jet sep¬ 
arately. The pc-jet power considered is the power estimated at 
the jet base. The blob radiation density is composed of densities 
radiated by synchrotron, SSC, second-order SSC, EIC from the 
BLR, and EIC from the pc-jet (t/; ^ -b t/; -t t/; + 


Table 4. Powers of the different components of the Ap Lib total 
jet expressed in log(P [erg.s^']) 


Power 

Blob 

Jet 

Total 

Radiation 

42.7 

41.7 

42.7 

Magnetic 

40.9 

41.2 

41.4 

Cold electrons 

43.5 

42.6 

43.5 

Non-thermal electrons 

43.9 

42.0 

43.9 

Protons 

46.8 

45.5 

46.8 


eicBLR eicj^' radiation density is the sum of energy 

densities radiated by synchrotron and SSC (t/'. = +t/'„_). 

The source is dominated by the kinetic power of particles, 
mainly in the cold protons of the blob, which represent 99.8% 
of the total power (see Table 4). The non-thermal electrons 
slighty dominate the cold electron population in the blob, while 
it is the opposite in the jet, but the two cases show almost 
an equipartition between cold and non-thermal electrons. The 
magnetic held appears far below the equipartition with the 
non-thermal electrons and the radiation in the blob, while this 
equipartition is almost achieved in the longer lifetime jet. 

Figure 7 adapted from Ghisellini et al. (2011) allows high¬ 
lighting the peculiarity of Ap Lib. The right part of the graph 
(Ld > 10"^^erg.s“') shows the very bright sources, mostly of the 
FSRQ type. In this part, Pr is proportional to the disk luminosity 
Ld. Assuming that Pr is also proportional to the jet power Ptotai, 
this means that for these sources the accretion efficiency, hence 
the accretion regime, remains constant. This is not the case any¬ 
more for disk luminosities below 10^^ erg.s ' where the radia¬ 
tion power instead evolves in VZ^. These two slopes of the ac¬ 
cretion regime correspond to two different accretion disks mod¬ 
els. On the one hand, the most luminous sources are described 
by standard disks, on the other hand, fainter sources are more 
consistent with disk models called advection-dominated accre¬ 
tion flow (ADAF), or radiative inefficient accretion flow (RIAF). 

The jet radiation power should normally place Ap Lib as BL 
Lac in the ADAF regime. However, although it is below this 
critical disk luminosity, Ap Lib seems to better follow the FSRQ 
type accretion regime, which is located in the prolongation of 
the grey FSRQ band, as could be an extremely faint FSRQ. One 
possibility could be that Ap Lib is an old FSRQ with an ac¬ 
cretion disk that became weaker and drifted toward the ADAF 
region (of low luminosities), but stayed meta-stable in the stan¬ 
dard disk regime. This idea is supported by all the signs high¬ 
lighted in our present study that are typical for an FRII type 
source such as the superluminal radio knots, the low aperture 
of the jet, and the large kinetic power dominance. In the follow¬ 
ing we also check whether the accretion efficiency is consistent 
with an FRII source. 

6.4. Accretion efficiency 

Since the kinetic protons power widely dominates all other com¬ 
ponents, the mass outflow rate can be expressed as 

Mou, = ^ - 6.0 X 1024^.5-'. (28) 

Taking in account that Mout < Min, with Mi„ the mass inflow 
rate (or mass accretion rate), allows placing an upper limit on the 
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Fig. 7. Blazar population in terms of disk luminosity Ld and jet 
radiation power Pr (figure in Ghisellini et al. (201 1)). The large 
red point represents the location of Ap Lib. The blue dotted line 
represents the separation between the two accretion regimes at 
Ld - 10"^^ erg.s The radiation power considered here for Ap 
Lib is the sum of the blob and jet radiation powers. 


accretion efficiency rj — Ldl{Mi„c^), 

i]< =93QxlQ-\ (29) 

^OUt^ 

This efficiency is lower than the one generally used for 
quasars of 77 ^ 0.1 (Raimundo et al. 2012). For Mom - Mm, the 
accretion efficiency is already at the limit between ADAF and 
standard disk regime of 77 ^ 0.01 proposed by Wu & Cao (2008). 


Ghisellini &. Tavecchio (2008) classified different blazar 
types in terms of rjcrit - LdjL^ddi which can be identified as 
Mm - Merit, with Merit - LEddlc^ the critical mass inflow rate. 
The Eddington luminosity of Ap Lib is 


Lsdci - 


AnGM,m„c 


O’ j 


M, , 
1.3 X 10^^— erg.s^', 

Mo 


(30) 


and its black hole mass is estimated as M = lO^ '^'^^^^Mo 
(Barth etal. 2003), which gives = L60+0.07X lO'^^erg.i^'. 
This corresponds to rjerit - 3.1 x 10“^, very close to the limit 
value of 3 X 10“^ between BL Lac and FSRQ proposed by 
Ghisellini & Tavecchio (2008). 


As shown in Sect. 6.3, the value Tjerit does not allow distin¬ 
guishing the various luminous FSRQs because the accretion ef¬ 
ficiency beyond Ld - 10"^^ erg.s^' seems constant. However, it 
makes a difference between BL Lacs and FSRQs, and also be¬ 
tween various BL Lac sources. 

These estimates of the accretion efficiency suggests that the 
source is in a peculiar accretion mode intermediate between 
ADAF and standard disk. Moreover, the radio features analysed 
in previous sections as well as the extended X-ray emission stud¬ 
ied by Kaufmann et al. (2013) show an intermediate blazar be¬ 
tween FR I and FR II types. In all aspects, Ap Lib corresponds 
to an intermediate case in the global blazar dichotomy scheme. 


which usually ties FRSQ with FR II and standard disk on the one 
hand and BL Lac with FR I and ADAF disk on the other hand. 

7. A specific type of LSP biazars? 

Simple one-zone SSC scenarios are often able to reproduce the 
SED and spectra of HBL sources reasonably well. In this type of 
blazar, the pc-jet, in which is embedded the VHE emitting blob, 
mainly contributes to the radiative output at low radio frequen¬ 
cies. At high energies, the contribution of this jet is weak or even 
negligible, for instance in the well-known HBL PKS 2155-304 
(Aharonian et al. 2005; Katarzynski et al. 2008). Conversely, the 
case of Ap Lib appears very uncommon in the sense that the pc- 
jet seems to dominate the SED in the radio and optical ranges 
and significantly contributes to the VHE part of the spectrum. 
Moreover, the strong radiation of the accretion disk places the 
source between the two populations of BL Lac and ESRQ. These 
properties make Ap Lib a rather unique object for the moment. 

To investigate whether Ap Lib is a special type of blazar or if 
it could represent a new still unrecognized family of sources, we 
have identified three main characteristics of its SED as follows: 

1) A broken radio-mm spectrum around 250 GHz. 

2) An X-ray positive slope. 

3) A high-energy bump very much wider than the syn¬ 
chrotron bump with a relatively flat Eermi spectrum at v < 10^^ 
GHz and a peak luminosity close to the peak synchrotron 
luminosity. 

We interpret the first characteristic as a superposition of two 
synchrotron bumps, one from a strong jet and one from a com¬ 
pact energetic blob embedded in it. The second characteristic 
is the signature of the inverse-Compton component dominance 
from soft X-rays to high energies. The third characteristic is ex¬ 
plained by a strong external inverse-Compton emission on the 
jet and disk radiation. 

Characteristics 2 and 3 are relatively common to de¬ 
scribe LBL spectra. However, characteristic 1 is very peculiar. 
Eollowing the description presented in Eig. (2) of Meyer et al. 
(2011), this change of slope at 250 GHz would be a key point 
for recognizing a source as a misaligned radio galaxy. But con¬ 
sidering that Ap Lib is a blazar with a very low angle with the 
line of sight (see Sect. 5.1), we explain this feature by an un¬ 
usual strong pc-jet emission and not by a strong misalignment. 
Therefore this radio-mm change of slope can be seen as a major 
point for characterizing the source. 

To see how common these three features are, we consid¬ 
ered the SED of 105 biazars obtained simultaneously by Planck, 
Swift, and Eermi reported by Giommi et al. (2012). Indeed, 
twenty of these 105 sources fulfil the three characteristics. All of 
them are classified as LSP BL Lac or ESRQ. Almost all sources 
with a spectral radio-mm slope change around 250 GHz have a 
positive X-ray slope and a wide Compton bump. Accordingly, 
we name here the sources that fulfil these three criteria BSRQs 
for broken-spectrum radio quasars, in which the emission of the 
radio VLBI base is relatively high and the blob synchrotron 
emission reaches maximum at very low-frequency. These two 
components significantly contribute and interact to produce the 
non-thermal emission from radio up to the VHE range. 

Thus Ap Lib apparently is not the only member of its type 
and can be seen as a prototype of VHE biazars, which could 
be seen as a radiatively extremely faint ESRQ, intermediate be¬ 
tween pure ESRQs and pure BL Lacs. In these types of specific 
VHE biazars, the whole SED is no longer dominated by a one- 
zone VHE blob as in the case of HBL sources, and additional 
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contributions from the disk and the jet become signihcant even at 
very high energies. However, a deeper analysis, which is beyond 
the scope of the present paper, is necessary to better characterize 
and identify these specihc types of LSP blazars and describe the 
transition between them and other blazar sources. 


8. Conclusion 

We conhrmed that a simple leptonic model dominated by the 
SSC emission of a single VHE compact zone (or blob) is unable 
to reproduce the SED of the blazar Ap Lib. We showed that addi¬ 
tional contributions related to the pc-jet and disk radiation appear 
to play a signihcant role, even at very high energies. This is con¬ 
sistent with the presence of an unusual strong jet emission, rein¬ 
forced by the observation at larger scale of an extended X-ray jet. 
To describe these different components self-consistently, we de¬ 
veloped a new code that includes the inverse-Compton emission 
of the blob particles on the jet radiation and on the disk radia¬ 
tion scattered by the BLR clouds, together with the absorption 
effects of the blob radiation by the jet. These effects naturally 
appear in blob-in-jet scenarios, but were previously found to be 
negligible for the modelling of HBL sources at very high ener¬ 
gies. However, specihc features of the SED of the LBL object 
Ap Lib can naturally be interpreted as their signatures and al¬ 
low constraining their physical parameters. On one hand, the two 
populations of relativistic particles, in the pc-jet and in the blob, 
provide a direct explanation for the spectral slope change de¬ 
tected in the radio-mm range, and enable reaching the observed 
VHE huxes thanks to the new inverse-Compton component on 
the jet radiation. On the other hand, the disk radiation and its 
subsequent scattering by the BLR clouds can simultaneously re¬ 
produce the spectral index seen in the UV band with UVOT and 
the apparent flatness of the Eermi spectrum. 

After identifying the pc-jet as an important component of the 
global radiative output of the Ap Lib nucleus, we constrained 
several of its characteristics from available VLBI data. The ob¬ 
served radio core properties constrain the aperture angle of the 
pc-jet and the angle with the line of sight, moreover, the veloc¬ 
ity and the size of radio knots observed in the jet constrain the 
knots expansion and the Doppler factor. We deduced a structure 
with an inner jet embedded in the pc-jet at small viewing an¬ 
gle to the observer. The simulated blob features are consistent 
with the radio knot observations, but the VHE blob is probably 
much closer to the core than the radio knots. The VHE blob be¬ 
haves like an emerging knot in expansion, which can later evolve 
into some VLBI knot and hnally contribute to form the radio jet. 
The source appears dominated by the kinetic power of particles, 
mainly by the cold protons of the blob. 

The classihcation of Ap Lib as a BL Lac seems not obvi¬ 
ous anymore. Its radio morphology, SED shape, accretion disk 
luminosity, powers, and X-ray properties correspond to an in¬ 
termediate case between ESRQs and BL Lacs. Although Ap 
Lib is a uncommon object among VHE blazars, the peculiar¬ 
ities found in its SED appear not unique. The spectral index 
change in the radio-mm range around 250 GHz is found in sev¬ 
eral other LSP BL Lacs or ESRQs that share other SED simi¬ 
larities with Ap Lib, namely an X-ray positive slope and a large 
Compton bump. We proposed to name this peculiar case BSRQ 
for broken-spectrum radio quasar, in reference to this unusual 
radio-mm index change in the SED and to distinguish them from 
the standard LBL category, which is composed of ER I type 
sources without such a strong pc-jet. An in-depth study of these 
peculiar BSRQs blazars is necessary to achieve a better under¬ 


standing of their mechanisms and the classihcation, but this is 
beyond the scope of this paper. 
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Note added in proof While this article was with the referee, the 
H.E.S.S. Collaboration et al. (2015) redefined the VHE spectrum of 
Ap Lib with more observation time. They reported a flux of (8.8 ± 
L5,„„, ± L8j^j)“'^ cm“^ s“' above 130 GeV , and a spectral index of 
r = 2.65 ± 0. 19 jm, ± 0.20jyj between May 10, 2010 and May 08, 
2011. Although it includes observations from 2011, this spectrum re¬ 
mains close to the one used here. It does not affect the capacity of the 
model developed here to reach the VHE range and therefore remains 
fully compatible with all results discussed above. 
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